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Abstract 


Viscous effects on transonic airfoil 
stability and response are investigated using an 
integral boundary layer model coupled to the 
inviscid XTRAN2L transonic small-disturbance 
code* Unsteady transonic airloads required for 
stability analyses are computed using a pulse 
transfer-function analysis including viscous 
effects. The pulse analysis provides unsteady 
aerodynamic forces for a wide range of reduced 
frequency in a single flowfield computation. 
Nonlinear time-marching aeroelastic solutions 
are presented which show the effects of 
viscosity on airfoil response behavior and 
flutter. Effects of amplitude on time-marching 
responses are demonstrated. A state-space aero- 
elastic model employing Pads approximants to 
describe the unsteady airloads is used to study 
the effects of viscosity on transonic airfoil 
stability. State-space dynamic pressure root- 
loci are in good overall agreement with time- 
marching damping and frequency estimates. 
Parallel sets of results with and without 
viscous effects reveal the effects of viscosity 
on transonic unsteady airloads and aeroelastic 
characteristics of airfoils. 


Nomenclature 


*h 



Cft 

C, h 

> 

C *6 


c "'h 

c m a 

c m 5 


= nondimensional distance from 
midchord to elastic axis 
= airfoil semi -chord 
^ airfoil chord 
= control surface moment 

coefficient about hinge axis due to 
control surface rotation 
= lift coefficient 
- lift coefficient due to plunge 
= lift coefficient due to pitch 
= lift coefficient due to control 
surface rotation 

= moment coefficient about pitching 
axis 

= moment coefficient about pitching 
axis due to plunge 
= moment coefficient about pitching 
axis due to pitch 

= moment coefficient about pitching 
axis due to control surface rotation 


k 

[K] 

m 


M 

[M] 

Q 


Qf 

Re 

s 

t 


pressure coefficient 
critical pressure coefficient 
nondimensional plunge displacement, 
positive downward from elastic axis 
dynamic plunge amplitude 
wb/U, reduced frequency 
structural stiffness matrix 
airfoil mass per unit span 
freestream Mach number 
structural mass matrix 

- (U/bw^j 2 , nondimensional 
dynamic pressure 
i (Up/bu^) 2 , nondimensional 

flutter dynamic pressure 

airfoil radius of gyration about 

elastic axis 

Uc/v, Reynolds number 

a + iu, Laplace transform variable 

time, sec. 



a m 

a o 

6 


U) F 

w h 


w a t , nondimensional time 
freestream velocity 
flutter speed 

distance aft of leading edge 
nondimensional distance from 
elastic axis to mass center 
airfoil angle of attack, positive 
leading edge up 
airfoil mean angle of attack 
airfoil dynamic pitch amplitude 
control surface deflection angle, 
positive trailing edge down 
control surface mean deflection 
angle 

control surface dynamic deflection 
amplitude 

m/Tipb 2 , airfoil mass ratio 
kinematic viscosity 
freestream air density 
Ut/b, nondimensional time 
angular frequency 
flutter frequency 
uncoupled natural frequency of 
plungi ng 

uncoupled natural frequency of 
pitching about elastic axis 


Introduction 


Research directed at better understanding 
of aerodynamic and aeroelastic phenomena at 
transonic speeds has increased greatly in recent 
years. These developments have been made 
possible by the advances made in computer power 
and numerical solution techniques. 1 Much effort 
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has been put Into the development of accurate 
yet cost efficient methods of calculating 
unsteady transonic airloads as well as the pre- 
diction of aeroelastlc characteristics such as 
flutter and divergence. 

Solutions for Invlscld unsteady transonic 
flowfields about oscillating airfoils have been 
made possible using computer codes such as 
LTRAN2. 2 While the Invlscld flowfields predict- 
ed by LTRAN2 give valid results for many cases, 
these solutions are not adequate when viscous 
effects are Important. Rlzzetta 3 Incorporated a 
viscous boundary layer model into the LTRAN2 
code using a viscous ramp model and Green's lag- 
entrainment equations, an Integral boundary 
layer model. A non-iterative Implicit technique 
was used to couple the quasi-steady viscous 
equations with the Invlscid potential flow equa- 
tions. GuruswaiTiy and Goorjian 14 applied this 
modified version of LTRAN2 to Investigate the 
effects of viscosity on transonic aerodynamic 
and aeroelastlc characteristics of oscillating 
airfoils. Two degree-of -freedom (plunge and 
pitch) typical section flutter speeds were 
presented for the conventional NACA 64A010 Ames 
airfoil (herein referred to as NACA 64A010A) and 
the supercritical MBB-A3 airfoil. As many as 
8000 time steps per cycle of oscillatory motion 
were required In the viscous lag-entrainment 
calculations to obtain reasonably converged 
results. Steady and unsteady viscous results 
agreed better with experiment than did compar- 
able invlscld solutions. 

Houwlnk 5 developed a procedure for an 
explicit coupling of the unsteady transonic flow 
and turbulent boundary layer computations, and 
Incorporated them Into the LTRAN2-NLR 6 code. 
Satisfactorily converged results were obtained 
In as few as 120 time steps per cycle using the 
resulting LTRANV 5 code. More recently, Howlett 7 
has modified the viscous coupling procedure of 
Rlzzetta for application to the XTRAN2L 8 general 
frequency transonic smal l -disturbance code. An 
integration of Green's lag-entrainment equations 
is performed from a specified transition point 
on the airfoil to the downstream boundary. The 
viscous ramp model Is discarded and boundary 
layer smoothing is used to reduce flow Instabi- 
lities. Iteration of an explicitly coupled vis- 
cous boundary layer solution with the Invlscld 
outer flow at each time step provides converged 
flowfleld solutions In as few as one or two 
Iterations. In general, converged oscillatory 
solutions Including viscous effects may be 
obtained In three cycles of motion with 360 time 
steps per cycle. Viscous solutions are obtained 
In computer times approximately twice that of 
the comparable Invlscld solution. Therefore, a 
more practical and affordable assessment of the 
effects of viscosity on transonic unsteady 
aerodynamic forces and aeroelastlc characteris- 
tics of oscillating airfoils Is now possible. 
The explicitly coupled vlscous-invlscld proce- 
dure of Ref. 7 Is used In the present study to 
calculate time-marching flutter solutions and 
transonic unsteady airloads required for aero- 
elastlc stability analyses. 

In the present study, airfoil stability 
analyses are performed using state-space aero- 
elastlc modeling such as that reported In Refs. 
9-12. Edwards, et al. 9 used a state-space model 


employing Padd approxlmants to model the 
unsteady airloads and demonstrated good agree- 
ment with a time-marching technique for a 
linearized case. Bland and Edwards 10 
demonstrated that such locally linear procedures 
may be used with airloads derived from a tran- 
sonic small -disturbance code. Batina and 
Yang 11 * 12 used a similar procedure to study 
the transonic aeroelastlc stability and response 
behavior of airfoils with active controls. 

The purpose of this paper Is to further 
Investigate the effects of viscosity on tran- 
sonic airfoil stability and response. The 
objectives of the study were: (1) to Investigate 
application of the pulse transfer-function 
analysis of Seidel, Bennett, and Whitlow 13 to 
treat airloads Including viscous effects; (2) to 
Investigate application of the nonlinear time- 
marching flutter solution procedure of Edwards, 
et al. 9 to transonic airfoil response Including 
viscous effects; (3)' to determine the effects of 
amplitude on Invlscld and viscous time-marching 
flutter solutions; (4) to assess the accuracy of 
state-space aeroelastlc modeling to predict air- 
foil stability Including viscous effects by com- 
parison with time-marching analyses; and (5) to 
apply state-space modeling to study the effects 
of viscosity on transonic airfoil stability. 
Results are presented for the NACA 64A010A, 
MBB-A3, and NACA 64A006 airfoils. Time-marching 
response and state-space stability analyses are 
performed for the NACA 64A010A and MBB-A3 
airfoils. Parallel sets of results are present- 
ed with and without viscous effects to determine 
the effects of viscosity on transonic airfoil 
stability and response. 


Computational Procedures 
XTRAN2L Transonic Code 

The finite-difference code XTRAN2L solves 
the transonic smal 1 -disturbance potential equa- 
tion and thus provides predictions of unsteady 
transonic flowfields about oscillating airfoils. 
The Invlscid algorithm, however, tends to over- 
predict the shock strength and locates the shock 
too far aft for cases where viscous effects are 
important. In order to predict both shock 
strength and location more accurately, viscous 
corrections are needed. The effects of 
viscosity are accounted for by an Integration of 
Green's lag-entrainment equations along the air- 
foil from a transition point selected at 10% 
chord to the downstream boundary. Iteration of 
the vlscous-invlscid solution may be performed, 
ensuring converged flowfleld results at each 
time step. A 


Pulse Transfer-Function Analysis 

Unsteady aerodynamic forces required for 
stability calculations are computed using the 
pulse transfer-function analysis available In 
XTRAN2L. This analysis Is based on the assump- 
tion that the unsteady forces are locally linear 
about the nonlinear transonic mean flow. In the 
pulse analysis, the airfoil Is given a small 
prescribed pulse In a given mode of motion and 
the aerodynamic transients are calculated. For 
pitch motion the pulse Is given by 
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elastic transients using complex exponential 
( 1 ) functions of the form 


where At is the nondlmenslonal time step. 
Similar expressions describe the pulses for 
plunge and control surface motions. The 
unsteady aerodynamic transfer-function Is then 
determined by dividing a fast Fourier transform 
(FFT) of the output force time history by the 
FFT of the Input pulse motion. The pulse analy- 
sis provides unsteady aerodynamic forces for a 
wide range of reduced frequency in a single 
flowfield computation. This is in contrast to 
multiple flowfield computations required for 
calculating oscillatory forces for discrete fre- 
quencies. Both invlscld and viscous pulse 
transfer-function analyses are performed. A 
typical Invlscld pulse analysis was computed 
with at = 5 it/32 and 1024 time steps while a 
typical viscous pulse analysis was computed with 
at = 2tt/27 and 2160 time steps. 

To assess the accuracy of the pulse 
transfer-function analysis Including viscous 
effects, unsteady aerodynamic forces were 
calculated for simple harmonic motion at several 
discrete values of reduced frequency k. Three 
cycles of oscillatory motion with 360 time steps 
per cycle and one viscous Iteration per time 
step were used for most cases. For k a 0.05, 
720 time steps per cycle and two Iterations per 
time step were required. 


Time-Marching Analysis 

The aeroelastlc system considered consists 
of plunge and pitch degrees-of -freedom. The 

equations of motion may be written In matrix 
form as 

{&*} < 2 > 

where the dot denotes differentiation with 
respect to nondimensional time t. Time- marching 
aeroelastlc solutions of Eq. (2) are obtained 
both with and without viscous effects. In the 
time-marching analysis, the equations of motion 
are coupled with the aerodynamic solution 
procedure of XTRAN2L for simultaneous 
time-integration. Equation (2) Is numerically 
Integrated in time using the modified state- 
transition matrix integrator of Edwards, et al. 9 
Details of the solution procedure may be found 
in Ref. 9. A typical Invlscld analysis was 

performed with At = 0.0004; a typical viscous 
solution was obtained with At * 0.0003. In 
general, several aeroelastic transients were 
calculated for a range of nondlmenslonal dynamic 
pressure Q. Values of Q were selected which 
resulted in subcrltlcal damped responses and 
supercritical diverging responses. The non- 
dimensional flutter dynamic pressure, Qp, was 
calculated by Interpolation and then confirmed 
by obtaining neutrally stable responses. 

Damping and frequency of the aeroelastlc 
modes are estimated from the transient response 
curves using the method of Bennett and 
Desmarals. These modal estimates are deter- 
mined by a least squares curve-fit of the aero- 
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Damping and frequency estimates from the time- 
marching analysis are plotted In the complex 
s-plane. 


Padd Model 


Aeroelastlc stability analyses are 
performed using a state-space aeroelastlc model, 
termed the Padd model, similar to that of 
Refs. 11 and 12. The locally linear Padd model 
provides a relatively Inexpensive determination 
of airfoil stability while retaining the non- 
linear properties of the mean flow. Transonic 
unsteady airloads required by this analysis are 
determined by Invlscld and viscous pulse 
analyses. The Padd model Is derived by assuming 
linear superposition of airloads due to airfoil 
plunge and pitch motions. The required airloads 
are approximated by curve-fitting the XTRAN2L 
unsteady aerodynamic forces with a Padd approxi- 
mating function. 11 The function may then be 
rewritten as a set of ordinary differential 
equations, which when coupled to the equations 
of motion, Eq. (2), and Laplace transformed 
leads to a linear first-order matrix equation 

l 1*1 - MU} (4) 

a 

where {z} contains the displacements, veloci- 
ties, and augmented states, and [A] Is a real 
matrix of constant elements. Equation (4) Is 
solved using linear eigenvalue solution techni- 
ques for specified values of Q. The resulting 
eigenvalues are plotted In a dynamic pressure 
"root-locus" type format and are compared with 
time-marching damping and frequency estimates In 
the complex s-plane. 


Results and Discussion 


Calculations were performed for the NACA 
64A010A, MBB-A3, and NACA 64A006 airfoils. 
Airfoil coordinates required for aerodynamic 
computations were taken from Ref. 15. Computa- 
tional conditions were selected to match the 
experimental conditions reported In 
Refs. 16-18. Transonic steady pressure 
distributions and unsteady aerodynamic forces 
were studied for all three airfoils for the 
computational conditions listed In Table 1. 
Time-marching response and Padd model stability 
analyses were performed for the NACA 64A010A and 
MBB-A3 airfoils. 

Aeroelastlc results are presented for two 
example sets of structural parameter values. 
Example 1 Is Case A of Isogal 19 which has bend- 
ing and torsion modes similar to those of a 
streamwise section near the tip of a sweptback 
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Table 1 Airfoils and computational conditions for 

transonic aerodynamic and aeroelastlc analyses. 


Airfoil 

Computational 

Condition 

M 

a m 

5 m 

Re 

NACA 64A010A 

1 

0.796 

- 0 . 21 ° 


12.56 x 106 

MBB-A3 

2 

0.765 

1.50° 


6.0 x 10 6 

MBB-A3 

3 

0.7557 

1.30° 



NACA 64A006 

4 

0.85 

0 ° 

0 ° 

2.41 x 106 


wing. Example 2 is the structural configuration 
used by Edwards, et al . 9 (also termed Example 2 
in Ref. 9) which has bending and torsion fre- 
quencies that are close together. Values for 
the structural parameters for the two examples 
are listed In Table 2. 


Transonic Steady Aerodynamic Results 

Steady flowflelds were computed for use as 
initial conditions for unsteady aerodynamic 
calculations. Steady pressure distributions for 
the lower surface of the NACA 64A010A airfoil 
are shown in Fig. 1 along with a plot of the 
airfoil contour. Both Invlscid and viscous 

XTRAN2L pressure distributions are presented and 
compared with the experimental data of Davis . 16 
Calculations were performed at the experimental 
conditions 16 of M = 0.796, = -0.21°, and 

Re - 12.56 x 10 6 which are herein termed Compu- 
tational Condition 1, as listed In Table 1. 
All three sets of results are In reasonable 
agreement. In the region of the shock, the 
viscous computations are in better agreement 
with the experimental data than the Invlscid 
computations. The calculated viscous shock Is 
slightly weaker In strength and Is located 
slightly upstream from the Invlscid shock loca- 
tion. Viscous effects are relatively mild for 
this case. 

Steady pressure distributions for the 
MBB-A3 airfoil are shown In Fig. 2 along with a 
plot of the airfoil contour. Computational 
results are obtained at two different sets of 
conditions for comparison with the experimental 


Table 2 Structural parameter values 
for aeroelastlc analyses 


Structural 

Parameter 

Example 

1 

2 

ah 

- 2.0 

-0.042 

*a 

1.8 

-0.036 


1.865 

1.368 

u 

60.0 

60.0 

“h 

100.0 

23.5 

“a 

100.0 

35.0 


data of Bucclantlnl, et al . 17 The first set 
corresponds to the uncorrected tunnel 
conditions 17 of M = 0.765, = 1.5°, and Re 
- 6.0 x 10 6 which are herein termed Computa- 
tional Condition 2, as listed in Table 1. The 
second set corresponds to the design condition 17 
of M = 0.7557 and a m = 1.3° which Is herein 
termed Computational Condition 3, as listed in 
Table 1. Computational Condition 3 allows 
Invlscid computations to match the experimental 
pressure data and hence viscous calculations are 
not presented at this condition. Both Invlscid 
and viscous solutions are presented at 
Computational Condition 2 . As shown in Fig. 2 , 
all four sets of results agree well along the 
lower surface of the airfoil except near the 
leading edge. For the upper surface, the 
XTRAN2L pressure distributions at Computational 
Condition 2 Indicate an overprediction of the 
shock strength and location, although the vis- 
cous computation Is In slightly better agreement 
with experimental data than the invlscid compu- 
tation. 1 Differences between the viscous results 
and experiment may be attributed to the fact 
that the boundary layer model does not account 
for strong Interaction between the shock wave 
and the boundary layer. Invlslcld calculations 
performed at Computational Condition 3, however, 
show much better agreement with experiment than 
either of the calculations performed at Computa- 
tional Condition 2. 

Similar comparisons of steady pressure 
distributions for the MBB-A3 airfoil at Computa- 
tional Condition 3 have been reported by 
Guruswany and Goorjian . 4 Similar to the Invls- 
cid XTRAN2L results at Computational Condition 3 
shown In Fig. 2, the viscous LTRAN2 results of 
Ref. 4 are In good agreement with the experimen- 
tal data of Ref. 17. Differences between 
XTRAN2L and LTRAN2 steady pressure distributions 
are attributed to the different grids and tran- 
sonic scaling used In the two codes. 

Steady pressure distributions for the lower 
surface of the NACA 64A006 airfoil are shown in 
Fig. 3 along with a plot of the airfoil contour. 
Both Invlscid and viscous XTRAN2L pressure dis- 
tributions are presented for comparison with the 
steady experimental data of Zwaan . 18 Computa- 
tions were performed at the experimental 
conditions 18 of M ** 0.85, = 0°, = 0°, 
and Re = 2.41 x 10 , which are herein termed 
Computational Condition 4, as listed In Table 
1. All three sets of results agree well except 
In the region of the shock. Here, the viscous 
computations are In better agreement with 
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Fig. 1 Steady pressure distributions for the 
lower surface of the NACA 64A01QA 
airfoil at M » 0.796 and - -0.21°. 



Fig. 2 Steady pressure distributions for the 
MBB-A3 airfoil at M - 0.765, * 1.5° 

and M * 0.7557, « 1.3°. 



Fig. 3 Steady pressure distributions for the 

lower surface of the NACA 64A006 airfoil 
at M - 0.85, * 0° and - OV 


experimental pressures than the invlscid 
computations. In the viscous solution, the 
shock Is weaker and Is located slightly upstream 
from the Invlscid calculation. Viscous effects 
are relatively mild for this case. 


Transonic Unsteady Aerodynamic Results 

Unsteady aerodynamic forces required for 
Padd model stability calculations were computed 
using the pulse transfer-function analysis. 
Representative results for the lift coefficient 
due to pitch, c* , are plotted In Fig. 4 as 
a function of reduced frequency k. These 
results were obtained using the pulse analysis 
with viscous effects for the NACA 64A01QA air- 
foil at Computational Condition 1 and a small 
pulse amplitude of ot 0 = 0.1°. To assess the 
accuracy of the pulse analysis, oscillatory 
calculations Including the effects of viscosity 
were performed for comparison. Results were 
obtained using an amplitude of a 0 = 0.1° at 
eight values of reduced frequency k = 0.05, 0.1, 
0.2, 0.4, 0.8, 1.2, 1.6, and 2.0. As shown In 
Fig. 4, • the pulse transfer-function Is In 
excellent agreement with the oscillatory air- 
loads. The excellent agreement between the 
two sets of results for this representative case 
clearly demonstrates the applicability of the 
pulse transfer-function analysis to Include vis- 
cous effects. 

To Investigate the effects of viscosity on 
transonic unsteady airloads, Invlscid and vis- 
cous pulse analyses are performed. A represent- 
ative comparison between pulse results for 
c i is given In Fig. 5 for the NACA 64A010A 
alFfoll (at Computational Condition 1). Both 
sets of results were computed using a pulse 
amplitude of 0.1°. The Invlscid and viscous 

pulse results show the same trends with respect 
to reduced frequency. Differences between the 
Invlscid and viscous c z results are largest 
for low values of k. a In the low k range, 
viscosity decreased the magnitude of both the 
real and imaginary parts. At high values of 
reduced frequency, the two sets of results are 
nearly the same. 

The effects of pulse amplitude on transonic 
unsteady airloads were Investigated by obtaining 
unsteady forces for successively Increased pulse 
amplitudes. Calculations were performed using 
both Invlscid and viscous pulse analyses. The 
Invlscid pulse amplitude computations (not shown 
here) Indicate that the unsteady forces are 
relatively Independent of amplitude. Viscous 
pulse amplitude computations showing effects of 
amplitude for Cg are presented In Fig. 6 
for the NACA 64A0?0A airfoil (at Computational 
Condition 1). Comparison between viscous pulse 
computations at a 0 = 0.1° with the computa- 
tions at a 0 = 1.0° shows some differences over 
the range of k plotted. With Increased pulse 
amplitude, the magnitudes of the real and Imagi- 
nary parts of c^^ are Increased. 

Unsteady experimental data from Davis, 16 
and the computational results of Houwlnk 20 and 
Guruswamy and Goorjian 4 are also plotted in 
Fig. 6 for further comparison. The experimental 
data of Ref. 16 was obtained using a harmonic 
pitch amplitude of a Q - i.o°. The unsteady 
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Fig* 4 Lift coefficient due to pitch about the 
quarter chord, c* , calculated using 
pulse and oscillatory analyses with 
viscous effects, for the NACA 64A010A 
airfoil at M * 0*796, % ■ -0.21°, 
and oq * 0.1°* 



Fig. 5 Comparison of Invlscld and viscous pulse 
results for the lift coefficient due to 
pitch about the quarter chord, c* » 
for the NACA 64A010A airfoil at 
M - 0.796, = -0.21°, and 

a 0 * 0.1°. 



Fig. 6 Comparison of pulse results with 
experimental data, LTRANV and 
LTRAN2-V1scous results for the lift 
coefficient due to pitch about the 
quarter chord, c* , for the NACA 
64A010A airfoil a? M * 0*796 and 

a m - -0.21°. 


forces of Ref. 20 were calculated with the 
LTRANV code using a 0 » 1.0°, M = 0.8, and 

a m = 0°. The unsteady forces of Ref. 4 were 

calculated with the viscous version of LTRAN2 at 
Computational Condition 1 using harmonic pitch 
amplitudes selected to match the experiment. As 
shown In Fig. 6, the viscous pulse computations 
performed using a 0 = 0.1° show good overall 
agreement with the experimental forces except 
for the Imaginary values of c^ at k = 0.1 

and k - 0.2. Comparison of the a viscous pulse 
computations with the LTRANV airloads shows good 
overall agreement. The imaginary part of the 

LTRANV results compares well with the Imaginary 
part of the viscous XTRAN2L results for pulse 
amplitude a 0 = 0.1°, whereas the real part of 

the LTRANV results are in better agreement with 
the real part of the viscous XTRAN2L results for 
pulse amplitude a 0 = 1.0°. Comparison of the 
viscous LTRAN2 results for c* with all of 

the other unsteady forces presented In Fig. 6 
shows good general agreement In the real part 
although the Imaginary part Is consistently 
overpredicted. 

Comparisons between invlscld and viscous 
pulse results for the unsteady aerodynamic 
coefficients due to control surface motion, 

c*-, c m , # and c hA» are presented In 
Figs. 7(a j, 7(b), and 7\c), respectively. These 
results were obtained for the NACA 64A006 air- 
foil with a trailing edge control surface of 25% 
chord (at Computational Condition 4). The 
Invlscld and viscous pulse results were 

computed using pulse amplitude 6 0 = 0.1°. The 
differences between invlscld and viscous results 
for c i and c m6 are generally largest 

for low values of k. In the low k range, the 

Inclusion of viscosity decreased the magnitude 
of both the' real and imaginary parts of c z - 
aind c m r . At higher reduced frequency, the 

two sets of results are nearly the same. In the 
pulse results for shown in Fig. 7(c), 

viscosity decreased the 0 value of both the real 
and Imaginary parts throughout the range of k 
plotted. Effects of pulse amplitude on 

Cfc 5 and c m(S are similar to those shown 
In Fig. 6. no amplitude effect was noted for 

c hs* 

Unsteady experimental data of Zwaan 18 and 
the computational results of Houwlnk 20 (calcula- 
ted using the LTRANV code) are also presented 
for further comparison (the harmonic control 
surface amplitude used In Refs. 18 and 20 was 
6p * 1«0°). The viscous XTRAN2L pulse solu- 
tion Is in better agreement with experiment than 
the Invlscld solution. Viscous pulse computa- 
tions show good agreement with LTRANV results 
for the real and Imaginary components of 
c A<5 and the Imaginary component of c^. 


Time-Marching Response Results 

Time-marching calculations were first 
performed using a small initial plunge displace- 
ment of h (0) = 0.001 to determine the value of Q 
which results In neutrally stable aeroelastlc 
transients. The effects of amplitude on flutter 
were then Investigated by obtaining time- 
marching responses for successively Increased 
Initial plunge displacement h(0). In these cal- 
culations, the dynamic pressure Q was set equal 
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(a) lift coefficient due to control 
surface motion. 

In viscid 



(b) moment coefficient due to control 
surface motion. 


Inviscid 



(c) hinge moment coefficient due to 
control surface motion. 

Fig. 7 Comparison of control surface pulse 
results with experimental data and 
LTRANV results for the NACA 64A006 
airfoil at M * 0.85, * 0°, and 

- 0 °; 


to the flutter dynamic pressure Qp previously 
determined using h (0) « 0.001. Only representa- 
tive responses for the NACA 64A010A airfoil (at 
Computational Condition 1) are shown In this 
section for the two example sets of structural 
parameter values listed In Table 2. 

Example 1. - Neutrally stable plunge and 
pitch time-marching response histories are 
presented in Fig. 8. The responses have been 
normalized by the Initial plunge displacement 
h (0) = 0.001. Inviscid and viscous flutter 
dynamic pressure values used to calculate the 
aeroelastlc transients are Qp - 0.60 and 

Qp = 0.97, respectively. The responses are 
bending dominated and are of constant amplitude 
after the higher frequency torsion mode tran- 
sients, visible in the first 2-3 cycles of 
motion, have damped out. As shown in Fig. 8, 
the viscous responses have the same characteris- 
tics as the inviscid responses Including nearly 
Identical amplitudes. Accurate modal curve-fits 
of the viscous response histories were obtained 
using the method of Ref. 14. Damping and fre- 
quency estimates of the aeroelastlc transients 
are plotted in the complex s-plane and are 
discussed in the following section. 

The effects of amplitude on the inviscid 
and viscous time-marching flutter solutions of 
Fig. 8 are shown In Fig. 9. Plunge and pitch 
degrees-of -freedom exhibit similar characteris- 
tics and hence only the pitch responses are 
shown. Furthermore, the pitch responses have 
been normalized by h(0) to allow for direct 
comparison between amplitude results. The pitch 
results for h(0) = 0.001 in Fig. 9 are identical 
to the pitch results of Fig. 8. As shown In the 
top part of Fig. 9, the inviscid pitch responses 
become slightly divergent when the Initial 
plunge displacement Is Increased by a factor of 
100. The amplitude effect is consistent with 
similar results reported in Ref. 9, where 
increased amplitude had a small destabilizing 
effect on the flutter responses of the NACA 
64A010A and MBB-A3 airfoils for the structural 
parameter values of Example 1. As shown in the 
bottom part of Fig. 9, the viscous pitch 
responses show a similar but larger diverging 
trend for h(0) = 0.05. At h(0) = 0.1, viscous 
responses diverged rapidly and led to program 
failure. The inviscid pitch responses of Fig. 9 
show a weak dependence on amplitude 
while the responses for the viscous case show a 
strong amplitude dependence. 

Example 2. - Neutrally stable time-response 
histories for the structural parameter values of 
Example 2 are presented In Fig. 10. These 
responses are of constant amplitude and 
oscillate about an exponentially decaying mean. 
In contrast with Example 1, the Inviscid and 
viscous flutter dynamic pressures are very 
close, Qp = 0.88 and Qp = 0.86, respective- 
ly. The inviscid and viscous plunge response 
histories have approximately the same amplitude 
and frequency although the viscous 
response oscillates about a slightly different 
mean than the inviscid response. Similar trends 
are visible in the pitch responses of Fig. 10 
(which are of smaller amplitude in comparison 
with the pitch responses of Example 1 (Fig. 8)). 
In contrast with the responses of Example 1, the 
higher frequency transient Is not present in the 
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Fig. 8 NACA 64A010A neutrally stable 

time-marching response histories of 
Example 1 at M - 0.796 and 
- - 0 . 21 °. 


Fig. 10 NACA 64A010A neutrally stable 

time-marching response histories of 
Example 2 at M ■ 0.796 and 
a m - -0.21°. 


Without viscous effects (Q f = 0. 60) 
h(0) = 0.001 


Without viscous effects (Qp = 0.88) 


-h(0) = 0.001 

hio)=ai 
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With viscous effects (Qp = 0. 97) 


With viscous effects (Qp = 0.86) 

h(0) = 0.001 

•2r- 

h(0) = 0.001 


Effects of amplitude on NACA 64A010A 
time-marching flutter solutions of 
Example 1 at M * 0.796 and 


Fig. 11 Effects of amplitude on NACA 64A010A 
time-marching flutter solutions of 
Example 2 at M « 0.796 and 




responses of Example 2 . The higher frequency 

mode is highly damped and thus has a negligible 
contribution to the total response. 

The effects of amplitude on the NACA 
64A010A time-marching flutter solutions (at 
Computational Condition 1) are shown in Fig. 
11. Inviscid pitch time-responses which are 
neutrally stable at h(0) = 0.001 remain neutral- 
ly stable at h(0) = 0.1. Viscous pitch 

responses are also neutrally stable when the 
initial plunge displacement is increased by a 
factor of 100. In contrast with the Example 1 
responses of Fig. 9, the Example 2 responses of 
Fig. 11 indicate that the flutter solution is 
relatively independent of amplitude. 


Pade Model Stability Results 

Pad6 model stability calculations were 
performed using transonic unsteady airloads 
determined by both inviscid and viscous pulse 
analyses. Pulse amplitudes were h 0 = 0.001 
and a 0 = 0.1°. Results for the NACA 64A010A 
and MBB-A3 airfoils are presented in this sec- 
tion for the two example sets of aeroelastic 
parameter values listed in Table 2. Comparisons 
between time-marchi ng and Padd model flutter 
solutions for Example 1 and Example 2 are 
presented in Tables 3 and 4, respectively. 

Example 1 . - Padd model dynamic pressure 
root- loci for the NACA 64A010A airfoil (at 
Computational Condition 1) are shown in Fig. 
12(a). Inviscid and viscous root-loci are 
plotted. With increasing dynamic pressure the 
torsion dominated mode moves directly left in 
the stable left-half plane while the bending 
dominated root becomes the flutter mode. The 
inclusion of viscous effects increased damping 
in the bending mode and slightly lowered the 
torsion mode frequency. Time-marching damping 
and frequency estimates are plotted in Fig. 
12(a) for Q = 0.4, 0.8, and 1.2. As shown in 
the figure, the Pads model root-loci are in good 
agreement with the time-marching modal estimates 
In both the inviscid and viscous cases. The 

value of dynamic pressure at flutter is given by 
the a = 0 crossing. The Padd model flutter 
dynamic pressure values are Qf = 0.62 and 
Qp = 1.00 for the Inviscid and viscous 


calculations, respectively. Thus, viscous 
effects increased Qp by approximately 61%. 
The Pads model Qp values are within 3% of the 
small amplitude time-marching values presented 
in the previous section. 

Pads model dynamic pressure root-loci for 
the MBB-A3 airfoil at Computational Condition 2 
are shown in Fig. 12(b). Both inviscid and 
viscous root-loci for the bending and torsion 
modes are plotted. Time-marching modal esti- 
mates are also plotted in Fig. 12(b) for 
Q = 0.2, 0.4, 0.6, and 0.8. Pads model results 
are in good overall agreement with time-marching 
damping and frequency values. Pads model 
dynamic pressures at flutter are Qp = 0.25 and 
Qp « 0.55 for the inviscid and viscous MBB-A3 
calculations, respectively. The inclusion of 

viscous effects increased the flutter dynamic 
pressure by approximately 124% in comparison 
with the inviscid calculation. Pads model Qp 
values differ from the small amplitude time- 
marching Qp values listed in Table 3 by 
7% and 22% for the inviscid and viscous cases, 
respectively. These differences may be attri- 
buted to the pulse amplitude effects on viscous 
transonic unsteady airloads similar to that 
shown in Fig. 6 and to the amplitude dependence 
of the flutter responses of Example 1 shown for 
the NACA 64A010A airfoil in Fig. 9. The larger 
differences between Padd model and time-marching 
Qp values in the viscous case may also be 
attributed to the more computationally sensitive 
nature of the viscous XTRAN2L calculations in 
contrast with the relatively routine inviscid 
calculations. 

The MBB-A3 root-loci of Fig. 12(b) are very 
similar to the NACA 64A010A root-loci of Fig. 
12(a). The inviscid flutter characteristics of 
these two airfoils were shown by Bland and 
Edwards 10 to be nearly identical when the steady 
shock strengths and locations were matched. The 
large differences in Qp between the two air- 
foils presented here may therefore be attributed 
to differences in steady shock strength and 
location as shown by comparison of Figs. 1 and 
2. Also, since the shock on the MBB-A3 airfoil 
(Computational Condition 2, Fig. 2) is stronger 
in comparison with the shock on the NACA 64A010A 
airfoil (Fig. 1), viscosity has more influence 


Table 3 Comparisons between inviscid and viscous time-marching and 
Padd model flutter solutions for Example 1. 


Airfoils and 

INVISCID 

VISCOUS 

Computational Conditions 

Time-marchi ng 

PadO Model 

Time-marching 

Padd 

Model 




Qp 

V“a 

Qp 

uv/u) 

F a 


u F /u) a 

NACA 64A010A 
M = 0.796, = -0.21° 

Re = 12.56 x 10 b 

0.60 

0.96 

0.62 

0.96 

0.97 

1.06 

1.00 

1.06 

MBB-A3 

M = 0.765, <»„, = 1.5° 
Re = 6.0 x 10 b 

0.27 

0.82 

0.25 

0.81 

0.46 

0.89 

0.55 

0.92 

MBB-A3 

M = 0.7557, a m = 1.3° 

0.65 

0.97 

0.74 

1.00 
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(a) NACA 64A01QA airfoil at 
M « 0.796 and am » -0.21°. 


(a) NACA 64A010A airfoil at 
M ■ 0.796 and 0 ^ * -0.21°. 


ODOfi Pad# model solution 
• ■+A Time-marching solution 



In viscid 

Viscous 



(b) M8B-A3 airfoil at M = 0.765 and 
“m - 1-5°. 


OQO^ Pad# model solution M = 0.765, o m a 15" 

• ■ + A Time-marching solution Mj5 0.7557, a *13* 



(c) MBB-A3 airfoil at M ■ 0.765, 
a m ■ 1.5° and M a 0.7557, 

«m * 1-3°. 

Fig. 12 Effects of viscosity on Pads model and 
time-marching (tynamlc pressure root- 
loci of Example 1; 


(b) MBB-A3 airfoil at M = 0.765 and 
a m = 1.5°. 


M = 0.765, a =15° 

m 

M= 0.7557, a m = 1.3° 



(c) MBB-A3 airfoil at M = 0.765, 
a m = 1.5° and M = 0.7557, 
a m = 1.3°. 

Fig. 13 Effects of viscosity on Pade model 
dynamic pressure root-loci of 
Example 2; 
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Table 4 Comparisons between inviscid and viscous time-marching and 
Padd model flutter solutions for Example 2. 


Airfoils and 

Computational Conditions 

INVISCID 

VISCOUS 

Time-marching 

Pad£ Model 

Time-marching 

Pad£ Model 

^F 

u F /u a 

■■ 

u F /u a 

mm 

“F /w o 

KB 

w F /w a 

NACA 64A010A 
M = 0.796, a m = -0.21° 
Re = 12.56 x 10 b 

0.88 

0.84 

0.97 

0.83 

0.86 

0.84 

0.92 

0.83 

MBB-A3 

M = 0.765, a m = 1.5° 
Re = 6.0 x 10 b 

0.61 

0,94 

0.58 

0.95 

0.86 

0.89 

0.96 

0.80 

M8B-A3 

M = 0.7557, a m = 1.3° 

0. 76 

0.83 

0.76 

0.83 
1 

i 

i 

! 

■ 


on the steady pressures and hence larger changes 
in Qp result. 

Inviscid dynamic pressure root-loci for the 
MBB-A3 airfoil at Computational Conditions 2 and 
3 are shown in Fig. 12{c). The Pade model 
flutter dynamic pressure values are Qp - 0.25 
at Computational Condition 2 and Qp = 0.74 at 
Computational Condition 3. The small decreases 
in Mach number and mean angle of attack in the 
inviscid calculations result in similar changes 
in the aeroelastic root-loci (Fig. 12(c)) as 
occurred with the inclusion of viscous effects 
at Computational Condition 2 (Fig. 12(b)). The 
three sets of computations performed for the 
MBB-A3 airfoil resulted in steady shock loca- 
tions of 63%, 61%, and 51% chord and monotoni- 
cally decreased shock strength as shown in the 
steady pressure distributions of Fig. 2. The 
corresponding Padg model flutter dynamic pres- 
sure values are Qp = 0.25, 0.55, and 0.74, 
respecti vely. The corresponding time-marching 
flutter dynamic pressure values are Qp = 0.27, 
0.46, and 0.65, respect i vely. With successively 
decreased shock strength and forward shock dis- 
placement, there is a monotonic increase in 
Qp. Pad£ model Qp values for the MBB-A3 
airfoil were typically nonconservative for 
Example 1 which may be attributed to the ampli- 
tude effects discussed earlier. 

Example 2 . - Inviscid and viscous Pade 

model dynamic pressure root-loci for the NACA 
64A010A airfoil (at Computational Condition 1) 
are shown in Fig. 13(a). In general, the inclu- 
sion of viscous effects increased damping in the 
lower frequency mode and decreased damping in 
the higher frequency mode. A switch in modal 
origin of the flutter mode occurs for this case 
with the addition of viscosity. The viscous 
computations give a flutter dynamic pressure 
value of Qp - 0.92 which is slightly less than 
the value of Qp = 0.97 given by the inviscid 
computations. The inviscid Pade model flutter 
dynamic pressure is within 10% of the time- 
marching flutter dynamic pressure, Qp = 0.88. 
Viscous Pads model and time-marching flutter 
dynamic pressures differ by only 7%, as listed 
in Table 4. The Pad6 model root-locus results 
were in qualitative agreement with the time- 
marching results. Damping and frequency 


estimates from the time-marching transients were 
not obtained, though, because of the closeness 
in frequency of the two aeroelastic modes. 

Pade model dynamic pressure root-loci for 
the MBB-A3 airfoil at Computational Condition 2 
are shown in Fig. 13(b). The inclusion of vis- 
cous effects increased damping and decreased 
frequency in the higher frequency mode. The 
lower frequency mode shows only small changes as 
a result of including viscosity. Dynamic pres- 
sures at flutter are Qp = 0.58 and 

Qp = 0.96 for the inviscid and viscous cases, 
respectively. Pade model Qp values differ 
from the small amplitude time-marching Qp 
values listed in Table 4 by 5% and 12% for the 
inviscid and viscous calculations, respectively. 
The inclusion of viscous effects resulted in a 
66% increase in flutter dynamic pressure, which 
is in contrast with the small changes in Qp 
found for the NACA 64A010A airfoil (Fig. 13(a)). 
The large increase in Qp between inviscid and 
viscous MBB-A3 cases at Computational Condition 
2 is due to the stronger steady shock in compar- 
ison with that of the NACA 64A010A airfoil (at 
Computational Condition 1). 

Inviscid Padg model dynamic pressure root- 
loci for the MBB-A3 airfoil at Computational 
Conditions 2 and 3 are presented in Fig. 13(c). 
The flutter dynamic pressure at Computational 
Condition 3 is Qp = 0.76 which is the same as 
the time-marching Qp value listed in Table 4. 
The small decreases in Mach number and mean 
angle of attack' in the inviscid calculations 
result in decreased damping In the lower 
frequency mode and increased damping in the 
higher frequency mode. A change in the modal 
origin of the flutter mode (Fig. 13(c)) also 
occurs which Is opposite to that for the NACA 
64A010A airfoil of Example 2 with the inclusion 
of viscous effects (Fig. 13(a)). Pads model 
flutter dynamic pressure values for Example 2 
were typically nonconservative which was similar 
to that found for Example 1. 


Concluding Remarks 

Viscous effects on transonic airfoil 
stability and response have been investigated 
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based on the use of the XTRAN2L transonic small - 
disturbance code. Aeroelastic results were 
presented for the NACA 64A010A and MBB-A3 air- 
foils. Aerodynamic calculations including vis- 
cous effects were performed using an 
integral boundary layer model coupled to the 
inviscid potential outer flow in a quasi-steady 
fashion. 

Transonic unsteady aerodynamic coefficients 
required for stability calculations were 
computed using a pulse transfer-function analy- 
sis. Excellent agreement was found between the 
pulse transfer-function and oscillatory airloads 
thus demonstrating the ability of the pulse 
analysis to include viscous effects. Inviscid 

and viscous pulse results showed the same trends 
with respect to reduced frequency although 

differences occur for low values of k. Viscous 
pulse computations showed better overall agree- 
ment with experimental data than the inviscid 
pulse computations. 

Nonlinear time-marching flutter solutions 
were obtained which showed the effects of vis- 
cosity and amplitude on airfoil response 

behavior and flutter. A two degree-of-freedom 
(plunge and pitch) aeroelastic system was con- 
sidered. Representative responses for the 

NACA 64A010A airfoil were presented for two 
different example sets of aeroelastic parameter 
values. In both examples, the viscous responses 
at flutter had the same characteristics as the 
inviscid responses including very similar ampli- 
tudes. In the first example, which had bending 
and torsion modes similar to those of a stream- 
wise section near the tip of a sweptback wing, 
large increases in flutter dynamic pressure Qp 
resulted with the inclusion of viscous effects. 
Also, the Inviscid flutter responses showed a 
weak dependence on amplitude while the viscous 
flutter responses showed a strong amplitude 

dependence. In the second example, which had 
bending and torsion frequencies that were close 
together, viscosity had little effect on flutter 
dynamic pressure. Also In the second example, 
the time-marching flutter responses were rela- 
tively independent of amplitude. 

Aeroelastic stability analyses were 

performed using a Pad£ state-space aeroelastic 
model. Dynamic pressure root-loci were present- 
ed for the NACA 64A010A and MBB-A3 airfoils for 
the two example sets of aeroelastic parameter 
values. Padg model root-loci for the 

first example considered, showed that viscous 
effects lowered the torsion mode frequency and 
increased damping in the bending mode thus 

delaying the onset of flutter. Changes In 
flutter dynamic pressure were correlated with 
changes in steady shock strength and location. 
Inclusion of viscous effects weakened the shock 
on the MBB-A3 airfoil more than that of the NACA 
64A010A airfoil and consequently caused larger 
increases in Qp. For the three sets of compu- 
tations performed for the MBB-A3 airfoil, 
decreased shock strength and forward shock dis- 
placement produced increased Qp values. In 
general, the Pad6 model root-loci were in good 
overall agreement with time-marching damping and 
frequency estimates in both the Inviscid and 
viscous cases. In the second example, viscosity 
caused an insignificant change in the NACA 
64A010A flutter dynamic pressure. For the 


MBB-A3 airfoil, a 66% increase In Qp resulted 
with the inclusion of viscous effects which was 
attributed to a stronger steady shock in compar- 
ison with the NACA 64A010A airfoil. Pade model 
flutter dynamic pressure values for both example 
sets of aeroelastic parameter values were typi- 
cally nonconservative. Differences between Pads 
model and time-marching values for Qp are 
attributed to pulse amplitude effects on viscous 
transonic unsteady airloads and to the amplitude 
dependence of the flutter responses In the first 
example considered. 
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